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Abstract 

We compute a model-independent correlation between the difference of neutron-neutron and 
proton-proton scattering lengths \ann — and the splitting in binding energies between Helium- 
3 and tritium nuclei. We use the effective field theory without explicit pions to show that this 
correlation relies only on the existence of large scattering lengths in the NN system. Our leading- 
order calculation, taken together with experimental values for binding energies and a^, yields 
ann = -22.9 ±4.1 fm. 
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I. INTRODUCTION 



In quantum chromodynamics (QCD) there are two effects that lead to violations of isospin 
symmetry. First, the electromagnetic interaction between quarks, and hence that between 
protons (p) and neutrons (n), does not respect the symmetry: is different from V^™ 

and V^^. Second, the difference between up and down quark masses {rriu^d) means that 
isospin symmetry is violated even in the absence of electromagnetic forces. The existence 
of TUd > rriu results in, for instance, different nucleon masses m„ > nip, which has profound 
consequences for nucleosynthesis in the early universe, and nuclear physics in general l| . A 
better understanding of isospin-symmetry breaking is therefore of deep interest to nuclear 
physics, and nuclear physicists. There has been much recent progress in this direction, with 
experiments at TRIUMF and lUCF exploring novel signatures of the violation of isospin 
symmetry, see Ref. for a review. 

The effect of isospin violation is significant in the nucleon-nucleon (NN) scattering lengths. 
That is because these scattering lengths are the result of fine tuning between the range and 
depth of the nuclear potential, and so small differences in either can lead to appreciable 
shifts in the scattering lengths. In a world of an isospin-symmetry-conserving interaction, 
both the neutron-neutron (nn) and the proton-proton (pp) channel are degenerate with the 
^5*0 neutron-proton {np, s) channel, because for two like fermions in a relative S wave, Fermi 
statistics allow for a spin-singlet, total spin J = state only. Hence, even relatively weak 
isospin- violating interactions could have a significant effect on the nn/pp system. 

The values of the proton-proton and neutron-proton scattering lengths, a^, Qnp^s are quite 
well established, see Table [H For charged particles, is observable as the leading-order 
(LO) parameter in a generalized effective-range expansion that includes the effects of the 
non-zero Coulomb repulsion in the asymptotic states. In contrast, pp phase shifts, and hence 
the pp scattering length, obtained exclusively from the strong part of the pp interaction 
retain a residual dependence on the specific model for this short-range force, i.e., they are 
dependent upon the renormalization scale and scheme sl-G]. 

Experiment reveals the differences between a^, a^p,,, and a^, and thus that charge 
independence and charge symmetry are both broken. Charge independence is associated 
with an arbitrary rotation in isospin space while charge symmetry is conserved if a rotation 
about 71 in isospin space leaves observables invariant (see, e.g. j^). Current data allows for a 
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relatively accurate extraction of np and pp scattering lengths, compared to the nn parameter 
where the uncertainty is about two orders of magnitude larger (see Tab. [I]). For a recent, 
thorough review of the status of experiments see Ref. Direct measurements are on the 
horizon js], but until now constraints on a„.„ come from final-state interactions. However, 
conflicting values for a„„ have resulted from attempts to follow this avenue in different 
few-nucleon reactions — indeed, conflicting values have resulted from different experiments 
investigating final-state interactions in n + d^n + n + p. 

To extract a„„ accurately from this three-body deuteron-breakup experiment, the outgo- 
ing particles should predominantly be in a state of very low neutron-neutron relative mo- 
mentum. Furthermore, the proton's effect on the detected neutron pair should be minimal. 
This condition is satisfied for a large separation of the scattered pair from the proton which 
remains at rest in the lab frame after the collision. This takes place in the nn quasi-free 
scattering (QFS) kinematics. The cross sections corresponding to this and the analogous 
np-QFS configuration constitute the experimental input for the subsequent extraction of 
ttnn- Currently, there are two data sets from which an identical theoretical method extracts 
conflicting values for ann- One setup records kinematical information of one neutron and 



the proton 



a 



10 



ing nucleons 11 



yielding a„„ = —16.1 ± 0.4 fm, while the other detects all three outgo- 



12| and produces a„„ = —18.7 ± 0.7. The theoretical model employs the 



CD-Bonn [isl NN potential and the charge-independent TM [l^ three-nucleon interaction. 
Recently, the sensitivity of theoretical predictions for the nn QFS cross section was investi- 
gated 15[. There it was shown that the nn-QFS cross section has a stronger dependence on 
rnn relative to changes in ann — after all, both experiments deal with small yet non-zero nn 
energies, and so r„„ would be expected to play some role. Ref. [l5| found values of a„„ and 
Vnn that plausibly fit both sets and would resolve the discrepancy, albeit at the expense of 
introducing appreciable charge-symmetry breaking in the nucleon-nucleon effective ranges. 

In view of these complications another way to "measure" ann would be of great interest. 
Here the effective field theory without explicit pions, EFT(/), is used to show that the 
difference of nn and pp scattering lengths, A(a) := a„„ — a^, is correlated with the trition- 
Helium-3 binding-energy difference, A (3) := B{t)—B{'^B.e). This correlation has been known 



for many yea rs within the context of models of the NN interaction. See, e.g. Refs. 16|-|20|, 



and Ref. 



2l|, which contains a review of work up until 1990. However, in these works 



various models, each of which produce a specific a„„, were used to compute A(3). The 
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possibility of mapping out the general relationship and using the result to constrain a„„ was 
not explored. Here we develop the correlation between A(a) and A(3) within a leading- 
order EFT{-/() calculation, which shows that this correlation stems solely from the existence 
of large scattering lengths in the NN system. We give an indication of how higher-order 
corrections can be expected to impact our result, and hence derive a constraint on a„„, 
using the value of given in Table HI 

In general, for systems where the scattering length a, is much larger than the range of the 
interaction R, an effective field theory based on the scale separation i? <^ |a| can be used 



to derive model-independent results 22N27j]. In nuclear physics this is EFT(/), and it is an 
expansion in R/a, with a given by the numbers in Table HI EFT(;<^) can be used to derive 
"universal" results that rely only on the existence of large scattering lengths. It has been 
used to compute triton {B{t)), Helium-4 {B{a)), and Helium-6 binding energies 28|-|30|. At 
leading order in the R/a expansion there are three parameters in the EFT that then yield 
predictions for all other observables in systems with A < 4: these can be taken to be the 
spin-singlet and spin-triplet NN scattering lengths, and the binding energy of the three- 
nucleon system. At next-to-leading order (NLO) in the R/a expansion the NN effective 
ranges enter the problem, with all other low-energy NN parameters only affecting answers 



beyond NLO 



31 



34j . This specifies the general EFT prescription of fitting a minimal set of 



LECs to observables in order to make predictions for all other observables correlated to the 
input set. 

In particular, EFT(;f) provides a map from an input set, e.g., the scattering length 
0'np,t, to a correlated set, whose elements, e.g., the deuteron binding energy B{d), are pre- 
dicted with known theoretical uncertainty. However, even if this set of correlations is well 
mapped out in the A-hodj system there appears to be no rigorous way to determine a 
priori whether A + 1-body observables will also be correlated with the A-hodj input quan- 
tities. For example, the triton binding energy B{t) is not correlated with the np singlet 
and triplet scattering lengths {anp...,, anp A (28[, while the binding energy of the a particle 
is correlated with {anp,s,cinp,t, B{t)} 29|]. For B(t), a strong sensitivity to short-distance 
structure — parameterized, for instance, by a momentum-space cutoff A — is found, while the 



dependence of B{a) on A is parametrical 



phenomenon, known as the Tjon line 
energy once that of the triton is known. 



y small once the value of B{t) is fixed. This latter 
35|, allows for a prediction of the a-particle binding 
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However, once it has been established that a higher-A observable is a member of the 
set of correlated quantities, we may use that observable to constrain properties of smaller 
subsystems. This prescription is used here, where we consider the binding energy difference 
A(3), which would be zero if isospin were an exact symmetry of nature. We therefore 
exploit another feature of the EFT, namely the absence of a qualitative hierarchy amongst 
low-energy observables in their role as input to fix the LECs. Thus, our interaction will take 
{ttnp^s, B{d),app, B{t), B(^B.e)} {A < 3) as input, and we will obtain the scattering length 
ann {A = 2) as output. Analogously, the Tjon line could be used to predict a range of B{t) 
values which are consistent with B{a). 

Our calculation of A(3) considers isospin violation from Coulomb interactions, and from 

n 

the difference in NN scattering lengths. Experimentally, A(3) is known to be 764 keV |57j . 
Note that we do not claim that our leading-order EFT(/) calculation of the individual tri- 
nucleon binding energies is this accurate, but we are examining a binding-energy difference 
that would be zero in the symmetry limit, and so a leading-order calculation of the difference 
already provides a useful constraint on |a„„, — a^|. In pursuing such a calculation we are, 
though, implicitly assuming that any isospin- violating component of the three-nucleon force 
in EFT(;f) enters only at sub-leading orders. We will present evidence that supports this 
assumption. 

The history of analyses of the impact of charge-independence breaking (CIB) and charge- 
symmetry breaking (CSB) on a„„ — and the trinucleon binding-energy splitting is a rich 
one (see e.g. 211] for a review of most of the investigations predating the advent of EFT 
methods in few-nucleon theory). Today, modern high-precision nucleon-nucleon force mod- 
els predict a trinucleon binding-energy difference, A(3), in good agreement with experiment: 
A(3,AV18+UIX) = 756(1) keV (l9| and A(3,AV18) = 762(9) keV |20|. Even though the in- 
dividual binding energies of the triton and Helium-3 receive significant contributions from 
a three-nucleon interaction (TNI) in these models, A (3) is driven by the difference in the 
nucleon-nucleon scattering lengths — at least once electromagnetic effects are properly ac- 
counted for. Isospin- violating TNIs were considered in the framework of chiral perturbation 



theory in Refs. 



36 



371 ]. and the leading-order isospin- violating TNI was found to contribute 



approximately 5 keV to A(3) [36|. 

A recent LO analysis of A(3) in EFT(^) included the Coulomb interaction nonperturba- 
tively. With a„„ and as input, a value of A(3) = 0.82 MeV was predicted using Faddeev 
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methods and the dibaryon formahsm 40|. In Ref. 4l|, the authors apphed EFT(/) at LO, 
NLO, and N^LO to predict proton-deuteron scattering- and bound-state observables. Their 
comprehensive analysis demonstrates the usefulness of the EFT prescription in low-energy 
pd scattering and in the Helium-3 bound state, since an order-by-order decrease of the theo- 



retical uncertainty is obtained. In Ref. j4l| too, a„„ and were used to fit the low-energy 



constants of the EFT. Neither of these EFT(/) analyses included isospin- violating TNI. 

The article continues with an introduction of EFT(;<') as the theory underlying the in- 
teraction potential. Next, a section on the numerical method of a refined version of the 
resonating group method (RGM), used to solve the few-body problem, precedes the pre- 
sentation of the results. The results section includes subsections discussing the uncertainty 
estimates due to suppressed higher-order long- and short-range interactions, and the limit- 
ing (hypothetical) case of a„„ — ^ o,pp- We then offer our conclusions, before assessing the 
numerical stability of the refined RGM in an appendix. 



II. INPUT: EFT(7f) WITH LEADING-ORDER COULOMB INTERACTIONS 
AND ISOSPIN VIOLATION 

The effective field theory without explicit pions (EFT(^)) at leading order in the expan- 
sion parameter Q/M is defined by the Langrangean 

C^""'^ = NUtdo + ^]N + CiiN^a2T2TiN) ■ {N^a2T2TiN) 

+C2{N^a2a,T2N) ■ (iVVaa.rsiV) + C^miN^ N){NW,N){N^nN) , (1) 

where the six-nucleon contact term renormalizes the S = 1/2 nucleon-deuteron channel 



(see, e.g. 28|]). The (iso)spin matrices (t)(T, with indices specifying the Cartesian com- 
ponent, project onto spin singlet and triplet with the respective low-energy constant, Ci^2- 
The two-nucleon amplitude derived from this Lagrange density matches the effective range 
expansion. The description is appropriate if the typical momentum exchange Q between 
interacting nucleons of mass is small relative to the high-energy scale M ^ mj^. For 
Q ^ JtItt, this theory is not applicable as is uses the neutron and proton Pauli spinors 
= (|p, s = 1/2) , |n, s = 1/2)) as degrees of freedom — nothing else. The interaction in 
Eq. ([1]) is charge independent and does not discriminate between neutron- neutron (nra). 
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proton-proton (pp), and proton- neutron (pn) pairs in the ^5*0 NN channel. For a compre- 
hensive analysis of systems including at least two charged protons at low energies, the effect 
of the electromagnetic force cannot be neglected. This is apparent in the measured difference 



between the pp and np ^Sq scattering lengths, = —7.8063 ± 0.0026 fm |5l| compared to 



anp,s = -23.748 ± 0.009 fm [5^, where Eq. (P yields a^^ = anp,s- 

Electromagnetic interactions are considered canonically by promoting the Lagrangean 
Eq. ([1]) to a local gauge theory, invariant under local U(l) transformations. Using Coulomb 
gauge V -A = 0, the contributions of the gauge fields to the pp amplitude can be split into 
a part which scales as a/Q"^, resulting from the part of the covariant derivative proportional 
to N'^cAqN ("Coulomb photons"), and others which are either suppressed by powers of 
Q/M or at least ("transverse photons"). Here, this estimate justifies the usage of the 
Coulomb potential resulting from the "exchange" of one Coulomb photon to account for the 
electromagnetic interaction at low energies. 

Without charge- independence-breaking (CIB) mechanisms stemming from a broken flavor 
SU(2) symmetry in the u-d quark sector, i.e., rriu 7^ rrid, the np and nn ^Sq channels would 
still be degenerate. To refine the analysis, the lowest-order contribution from this asymmetry 
is included. The dominating terms are expected to be of lowest-mass dimension while a 
dependence on the direction of the isovector is now admissible in order to distinguish nn, pp. 



and np vertices. In a hierarchy of isospin- violating interactions |42|, |55| , those contact terms 
which are expected to scale as eQ^ with e = ~ | should be subleading compared to 

the Coulomb potential. We include both the Coulomb potential, and these isospin-violating 
short-range operators, in our calculation, allowing us to formulate a model-independent 
assessment of which a„„ values are consistent with the experimental tri-nucleon binding 
energy splitting A(3), the singlet (triplet) np S-wave scattering lengths anp,s{t), and a^. 

The subsequent analysis is therefore based on the Lagrangean ([1]), including pieces ob- 
tained by the minimal substitution 9o — )■ 9o + ieAo, combined with explicit CSB terms: 

^(CSB) _ ^nn (^^tjy) (n^ a2 (l + iraTi) N) + Cf (N^N) (iVVs (1 - iT2n) N) . (2) 
The tree-level diagrams corresponding to these interactions define the isospin-violating part 
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of the potential 

V^(CSB) 



E 

Kj 



4|f(i) - f(j)| 

r3(z)) (1 - T,ij)) U{r-, 

A3 



+CT (1 



in coordinate representation with fhirij; .— y^:^^ 
conserving piece of the potential, which stems from Eq. 

A 



(l + r3(z))(l + r3(j)) 
1 



J 4 



(1 - ■ aU)) 



(3) 



e 4 



is: 



Meanwhile, the isospin- 



V^(ci) 



^(l-a(^)-a(j))Ci + ^(3 + a(0-a(j))C2 



(4) 



i<j<k 

cyclic 



These operators yC-^^^), V^'^^\ are used in the Schrodinger equation and originate from 
the momentum-independent vertices via a Gaussian regulator of the non-separable form 



yp-p 



7A2 



which Fourier transforms into a Gaussian depending on the relative 



coordinate. From this non-relativistic equation of motion, variational approximations to the 
bound and scattering states of the two- and three-nucleon systems are obtained. In solving 
the equation of motion, the potential is iterated. 
The iteration of the unregulated 



np total isospin T = amplitude 



38 



— oo) interaction ~ C2 in the ^Si channel yields a 
d^q 1 



1 



(27r)3E 



le 



(5) 



where p = \/m^E is the relative momentum of the np pair in the cm. frame and e is 
a positive infinitesimal. C2 can then be chosen in order to obtain the real bound state 
(deuteron) in that channel: 

T^=°--(p) = — (6) 
m^'^ + ip 

where we have chosen to ignore terms which are suppressed by p^/A in the denominator, 
although these pieces will be present in the denominator for any finite A, and will thus pro- 
duce a non-zero, positive 391], and cutoff-dependent effective range. The binding momentum 



is denoted here as 7 = ^ymj^B{d). 
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Similarly, the iteration of either Ci or Ci + C^"" produces the virtual bound states in the 
nn and np systems: 



1 



+ ip 



(7) 



under the same conditions as in Eq. ([6]). The isospin- violating terms account for the dif- 
ference in the scattering length of nn, relative to np. Without the C^" counterterm, the 
spin-singlet nn and np channel would be degenerate: a„„ = ttnp^s- 

The case in the proton channel is different, due to the presence of the Coulomb interaction 
there. The computation of pp scattering was carried out to LO in EFT(/) for S-wave NN 



scattering in Re: 
final result is 44 1 



431]. Our presentation of pp scattering rests on that treatment. The 



T — Tnc + Tcouh (8) 

with Tcoui the amplitude for scattering due to the Coulomb potential alone. T/vc encodes 
the purely strong scattering and the Coulomb-nuclear interference: 

1 



Tnc = C exp{2iao{ri))- 



- Mp) 



(9) 



Here is the Sommerfeld factor: 



2717] 



(10) 



and the Coulomb parameter rj := with (Ji{ri) = argr(/ + 1 + if]) (F is the Euler gamma 
function). In Eq. (j9]) Jo{p) is the Coulomb- modified bubble, depicted in Fig. [H If computed 
in PDS [231] at a renormalization scale /i, its finite part is: 



in 



\am^J 2 



47r ' 



(11) 



once divergences in D = 4 and D = ?> have been dropped. In Eq. ( ITTj) . the function 

H{'q) = i){ir]) + ln(2?7). 



(12) 

and 



with ip the derivative of the Euler Gamma function, Euler's constant Ce = 0.5772(. 
a = c^/Att (see also j46l-l48|). 

The denominator in Eq. (Q can be matched to the modified effective-range expansion: 

in exp(2i(Jo) 



-7VC 



(13) 



FIG. 1: A diagrammatic definition of the Coulomb bubble Jq. The solid lines are nucleons, and the 
wavy line represents a single insertion of the Coulomb potential. The shaded blob is the t-matrix 
found by iterating the Coulomb potential alone, Tcoui- 



This expression for T encodes the fact that the spherical Bessel functions, the asymptotic 
solutions for the np/nn system, are now replaced by Coulomb functions 4^. When this 
matching is performed at leading order (ro = 0) we find that 



ami f fJ.y/7r\ ^ ^ 3 



am^ I I 



(14) 

An ^ ' 



To produce the experimentally observed the combination C\ + C^^ is fitted to the value 
listed in Table HI Strictly speaking, the listed there is extracted from data using an 
effective-range expansion that considers a refined version of the Coulomb potential and 
vacuum-polarization effects 5l|. However, here we have a long-range part of the potential 
consisting only of the Coulomb interaction. The uncertainty due to higher-order long-range 
interaction effects, such as vacuum polarization, is included in the overall theoretical error 
estimate. 

Regardless of this detail though, we see that can be obtained from the amplitude 
in a manner that is independent of the renormalization scale In contrast, attempting to 
"switch off' Coulomb interactions and compute the effect obtained solely from the strong 
potential, C\ +C^^, yields a result that, within PDS, is dependent upon the renormalization 
scale yU. C\ + C^^ contains a In /i piece, and there is no corresponding In /i piece of the 
Coulomb- less loop function to cancel that. Our calculation is not done using PDS, but 
instead with a Gaussian regulator, yielding approximately a In A dependence of the C^g that 
reproduces (the In A dependence is seen explicitly with a sharp momentum cutoff {4]). 
But, in any case, this cutoff /renormalization-scale dependence renders quoting a "strong pp 
scattering length" a questionable exercise — at least in the absence of an agreed upon 
choice for the regularization and renormalization scheme and scale. Therefore, in what 
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follows we quote all results in terms of a^, which is a physical observable, and as such is 
independent of these choices. 
In practice, the equation 

^ + + ^ (^1-2, C3N1) j 1^) = E\ij) (15) 

is solved for two-body scattering (i? > , A = 2) and three-body bound states (i? < , A = 
3, see ch. IIIII below). In total, five LECs are fitted to low-energy data according to Tabled 
Again, it is the sum Ci + that controls the nn channel, not just C^"! 

Isospin-violating interactions at the quark level imply furthermore a mass difference be- 
tween the neutron and the proton. This consequence is a higher-order effect in our counting, 
as it is not enhanced by the fine-tuning in the NN system. Hence, it is not accounted for 
in the calculation below, where throughout m„ = rrip =: = ^(938. 211 + 939.505) MeV is 
used. In sec. lIVBt the contribution of this mass difference to A(3) is included in the error 
assessment. 

III. THE RESONATING-GROUP METHOD 

For the solution of the two- and three-body problem the variational Resonating Group 
Method (RGM) is employed. In the following, the ROM is introduced by the specific example 
of the variational space used for the triton calculation of this work. 

The ground state is expanded in two different variational bases which differ in the angular- 
momentum coupling scheme between the spin- and coordinate-space components of the wave 
function. First, the ansatz for the three-body state in the "LSJ-scheme" (total angular 
momentum J and parity vr), used here in Jacobi-coordinate space {pi,p2), reads 

1-0(^1, P2), J = 1/2, vr = +)^^ = fi{pi, p2)\d - n) +/2(pi,P2)|d - n) + h{pu P2)\{nn) - p) , 

(16) 

namely, a linear combination of all possible two-fragment substructures within the triton: 
deuteron (d), spin-singlet deuteron (d), and dineutron (nn). Thus, the coupling schemes of 
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the spin (S) and isospin (T) angular momenta are fixed, e.g., 



\d — n) 



|n, TUs^) ® |p, rris^) (g) |n, m^g) 



1 1 iO 



(17) 



The square brackets are shorthand for a basis where the irreducible representations are 
labeled by quantum numbers corresponding to {S+S')"^ for two general tensor operators with 
spherical components m(m'), [5"* ® 5""^^^ = ^„ ,„,(5m5W| JM)5'^5"^' (see e.g. |5o|). 
The coordinate-space wave functions fi are expanded in a Gaussian basis of dimension Df. 

Di 

U = Y^ Q, exp(-7rpi^ - 7?pt) ■ [3^0 (Pl) ® yo {P2t=' . (18) 



Here the two solid harmonics 



501 ] are both chosen to correspond to relative angular mo- 



28, 



411] with LO EFT(^), the dynamics in 



mentum zero. As in the Faddeev approach of 
the variable pi are affected by S-wave NN interactions only. The pi coordinate parameter- 
izes the two-nucleon fragment, for instance, the relative coordinate between a proton and 
a neutron in a deuteron-neutron configuration. As the lowest-lying bound state will occur 
for an angular momentum on the second coordinate, e.g., between the center of mass of a 
deuteron cluster and the remaining neutron, which is zero as well, the total orbital angular 
momentum is also zero (L = 0). The impact of higher angular momenta on the solution for 
the 3N bound state is assessed in Appendix Rl 

The relative size of the variational parameters Cij then determines the overlap of the cor- 
responding configuration with the triton, and hence the significance of a certain combination 
of width parameters '~fi2 ^^e grouping i. The triton ground state, for example, is found 
to have the largest overlap with the d-n grouping. To this configuration, parameters 71 
that resemble the spatial extent of the deuteron and 72 that place the neutron at a farther 
distance, contribute most. For the other two groupings, 2 = 2, 3, the two-body fragment is 
unbound. In consequence, basis vectors with smaller widths 71, i.e., broader spatial extent, 
become important. 

In the basis f|T6l) . the individual spin- and orbital angular momenta are coupled separately 
to their total S and L before they combine to total J. This coupling scheme is different 



from the one used in 



4l| . The RGM ansatz used for two-fragment scattering resembles that 
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TABLE I: Data input for the LECs which specify the interaction Eq. ([3]). The values were fitted 
nonperturbatively with the RGM and include the Coulomb potential's contribution. Four sets of 
LECs corresponding to A = 400, 800, 1600 MeV were adopted. The two sets for A = 400 MeV 
correspond to different input data. 



channel 



(np) ^So 
(np) 



{pp) 'So 
{npp) "^Si 
{nnp) 



low-energy constants (LEG) 



Ci + Cf 
Ci + Cg^, Ci , C2 , C3N1 

CI /^nn 
1 + 0_5 , Oi, O2, O3NI 



constraining observable 



anp,s = -23.748 ± 0.009 fm [52] 
a„p^t = 5.4194 ±0.002 fm 
Bid) = 2.224575 ± 0.000009 MeV [53] 
ag, = -7.8063 ± 0.0026 fm [51] 
S(3Hc) = 7.718109 ± 0.000010 MeV 
B{t) = 8.481855 ± 0.000013 McV 



calculation: 

|^/'(pi,p2),^= l/2,7r = +) 



ch 



/l 



\d) ® \n) 



-I J 



|d) ® \n) 



(19) 



with /j^j^^j = Xl^i Cije~^'^^'^yi^^^ (P2) and a two-fragment wave function determined within the 
same framework. Now, the deuteron's J = 1 is combined with the neutron to a channel spin 
Sch which couples to the relative orbital angular momentum ^='^ to the total J. With 
the two-fragment wave function fixed, the variational parameters Cy determine the relative 
importance of the components and the expansion of the coordinate-space wave function of 
the third particle relative to the center of mass of the pair. This resembles the usage of 



the various dibaryon propagators A^^l^^^^ for the deuteron (^^o) channel (see 28|, |4l|) in 



the Faddeev approach, and the corresponding projection to deuteron-neutron (for example) 
relative angular momentum zero. The expansion parameters Cij of the triton in both the 
LS and the channel basis are obtained by solving the generalized Eigenvalueproblem for a 
Hamiltonian incorporating the potential of Eq. ([3]), i.e. Eq. f|T5|l . 



IV. RESULTS 



Here, to obtain a result for ann different from the values of the other two-nucleon scattering 
lengths, the triton binding energy B{t) is used as input for C™. The triton constitutes 
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the next larger system where an interaction between two neutrons is observable. If there 
are no additional three-nucleon operators which affect low-energy modes and break isospin 
symmetry then B{t) and correlated, and therefore A (3) is correlated with A (a). An 

interaction given by Eq. with B(t) as input for C"'" will then predict a„„ within our 
uncertainty bounds. 

The results presented in Fig. |2] are RGM solutions to the Schrodinger equation with an 
interaction given in Eq. ([3]) and LECs constrained by the conditions defined in Table [B The 
correlation as shown there is obtained by a variation of C"". The predictions for a„„ of 
four interactions, each containing different short- distance physics (see below), inferred from 
the intersection of the computed correlation line with the experimental value for A(3) : = 
B{t) - 5(3He), suggest 

lim ann{A) ~ -20.8 fm. (20) 

A— >CXD 

An assessment of whether the chosen set of input quantities and leading-order calculation 
are such that this prediction discriminates between the two conflicting measurements of ann'- 



-18.7±0.7fm 
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12| and ann = -16.1±0.4 fm Q,Q occupies much of the rest of this 
section. If both data points are consistent with the LO prediction fl20l) within the uncertainty 
of that calculation, a higher-order analysis is needed before any conclusion can be drawn. 
That theoretical uncertainty results from two expansions and their respective truncation. 
The error from expanding the wave function in a finite-dimensional variational space (see 
above) is discussed in Appendix lAl and is found to be ±1.5 keV there. In Subsections IIV Al 
and IIVBI the effect of omitted higher-order interactions in the expansion of the Lagrangean 
shall be analyzed. The additional contributions result from the broken isospin symmetry 
of the underlying theory in both the strong and the electromagnetic sector. The omitted 
effects in the EFT are divided here into short-range parts, which all occur in the strong part 
of the charge-symmetry breaking Hamiltonian, and long-range effects, the most important 
of which are electromagnetic interactions beyond Coulomb repulsion and the impact of the 
neutron-proton mass difference on A(3). 

A. Theoretical uncertainty I (short-range) 

In this subsection the effect of omitted higher-order NN operators in the Lagrangean 
which correspond to strong interactions will be assessed. These operators must violate 
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FIG. 2: Correlation between the trinucleon-binding-energy difference B{t) — B(^]ie) and the nn and 
pp scattering lengths at LO in EFT(7f). Correlation lines correspond to regulator values 400 MeV, 
800 MeV, and 1.6 GeV. For A = 700 MeV, the correlations overlap with the displayed lines and 
are not shown. 

charge symmetry, or they will not contribute to A(3). They become relevant for modes with 
momenta of order m.^. In this bound-state calculation these high-energy modes can only 
affect observables through loops. Their effects are assessed here in two ways. First, by a 
change in the regulator parameter A; specifically, A = 400, 800, 1600 MeV and thereby three 
different short-range potentials are used. Second, by a change from anp,t to B{d) as input for 
6*2, which demonstrates how imposing renormalization conditions at different characteristic 
momenta changes the output predictions. 

The LEG values corresponding to these four different cases (A = 400, 800, 1600 MeV, 
anp,t instead of B{d) as input) are different because each set represents a different model 
for the short-distance part of the interaction. The substitution anp,t B{d) was made for 
A = 700 MeV only, and did not result in a significant shift of the correlation line. All 
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four sets of LECs reproduce, by construction, the same long-distance behavior, as given 
by the specific input data, within higher-order uncertainty. The difference in predictions 
derived from the four potentials is thus the uncertainty due to permitted variations in the 
short-distance physics. This sets a lower bound for the theoretical uncertainty. 

Our assumption that three-nucleon forces do not contribute to B{t) — B(^B.e) at leading 
order is supported by the absence of any significant cutoff variation in our result for this 
observable. This cutoff dependence is quantified by the width of the correlation band shown 
in Fig. [21 Naive dimensional analysis suggests that the three-nucleon contribution to this 
binding-energy difference should be less than that of CSB two-nucleon operators, and the 
results found here do nothing to contradict that view. Similar conclusions have recently been 



reached by other authors 



4l| . Of course, naive dimensional analysis is strongly violated in 



EFT(/) in the isospin-symmetric sector, with the TNI being of leading order there. Thus a 
power counting for isospin- violating operators in EFT^jt), e.g. along the lines laid out for 
sub-leading isospin-symmetric forces in Ref. 49|], remains an interesting open problem. 

The following brief discussion explains the choice of A = 1600 MeV as an upper bound 
for the cutoff. A is not increased further because of the corresponding increase in the 
unrenormalized three-nucleon binding energy. A diverging three-body ground-state energy 
is the result of a two-body interaction whose range is decreased and strength increased to 
fit cinp,s and anp,t (Thomas effect). The discrepancy in the binding energy with and without 
a three-nucleon force — introduced to properly renormalize one three-nucleon bound-state 
energy to that of the triton — increases in consequence, with more bound states entering the 
spectrum at specific cutoff values. In Fig. [3] the relevant part of the three-nucleon spectrum 
that is obtained without a TNI is shown as a function of the regulator. The NN interaction 
fits Ci.2 to the singlet and triplet np scattering lengths for each A, and a new "Efimov" 
trimer is found to enter around A = 750 MeV. The three-body force has then either to be 
strong enough to lift a very deeply bound three-body state (solid black line in Fig. [3j) to the 
experimental triton energy while unbinding the other states, or it must pull the shallowest 
of those states (solid red line in Fig. [3]) to the triton level. In either case the variational 
basis must be refined, either to expand states of considerably larger binding energy, or to 
treat excited states accurately. Convergence for either case requires a larger RGM space 
compared to the one where only one three-body state is bound, and bound with an energy 
already of the same order of magnitude as that of the triton. 
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FIG. 3: The two lowest eigenvalues of the nnp system of a Hamiltonian with the two-body inter- 
action fitted to anp,s and anp,t without a three-body term, as a function of the regulator width 



For the three representative cutoff values used to obtain Fig. [21 the three-body parameter 
was always adjusted to elevate the ground state to the triton level, because this repulsive 
force pushes all the shallow states that are not typically resolved by the RGM basis above 
the d — n threshold. (In the qualitative analysis of the previous paragraph, whose results were 
shown in Figl3l a small space was chosen that expands only one of the shallow states, EV2 
in Fig. [31) The long-distance properties of the three-body state are correctly reproduced by 
this basis, as witnessed by the fact that the d — n threshold approaches its predicted value 
of B{d) = 1 / {a^jj ^m^) ^ 1.41 MeV for A — )■ 00. (The value here is the correct one, given 
that anp,t is taken as input. Slightly different thresholds will be found at any finite A, since 
in that case higher-order terms in the effective-range expansion are not zero.) 

Bearing this limitation in mind, the graph in Fig. [2] allows the theoretical uncertainty at 
this order to be estimated as the maximal a„„ difference at the experimental binding-energy 
splitting: 



Accordingly, we identify the mean value of all regulator-dependent a„„ values with the 
EFT(^) prediction. We stress that this choice is somewhat arbitrary and does not coincide 



A. 




(21) 
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with the choice made in 



30| for the average charge radius of the triton. For that radius, the 



EFT(/) prediction was taken to be an extremal value calculated for the smallest employed 
cutoff. That choice was made because another LO calculation predicting the triton charge 



radius 



541 ] showed that the variation of cutoff and input data in 30| was not reflecting 



the total LO uncertainty. In fact, both here and in |30|] we take the mean of all available 
EFT(^) predictions for the observable of interest (here, ann and in 30|] the triton charge 
radius), regardless of the numerical method used to obtain them. The LO uncertainty is 
then given as half the difference between the minimum and maximum value predicted for 
that observable. This seems a sensible general prescription for defining central values, and 
uncertainties due to short-distance physics, in EFT(;f) calculations. 

Here, only the a„„ = —18.7 fm data point lies within this uncertainty range. This 
result, however, does not yet allow for the conclusion that the smaller datum is inconsistent 
with the input data, {a„p _B((i), 5('^He), a^}. Only if the uncertainty due to long-range 
contributions does not add to ^^hort an amount that would increase the total uncertainty to 
eventually include the second data point, can such a discrimination be made. 



B. Theoretical uncertainty II (long-range) 

As a cutoff variation assesses only the dependence on short-distance structure, an estimate 
of the theoretical uncertainty is incomplete without considering effects sensitive to low- 
momentum modes. Above it was argued why the long-distance part of the interaction can 
be approximated by the Coulomb potential. It is shown here that uncertainties from thereby 
omitted higher-order interactions, which are not assessed by the A variation considered in 
the previous subsection, are sufficiently small to render only one experimental data point 
for ann consistent with the input observables. 

For this purpose, A(a) is treated as a function of A(3), i.e., an interpolation of the 
dependence whose graph is shown in Fig. [2] is inverted. Furthermore, the assumption that 
higher-order interactions will yield a correlation line which might be shifted but is identical 
in shape compared to the ones shown is made. Then, the error in A (a) introduced by the 
suppressed terms who contribute to A (3) a correction Ac (3) can be approximated via 



gA(a) (A(3)) 
o.ong- a(A(3)) 
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■Ac(3) . (22) 

A(3)=0.756 MeV 



The two long-distance effects we must consider are: 

1. The difference in the nucleon kinetic energy due to non-equal neutron and proton 
masses. Parametrically we estimate the size of this effect to be: 

' )^(m„-m,)(-^). (23) 



2mAr 2m„ " ^ '^rn'j^ 

Taking {^) ~ B{t) we find: 



Bit) 

Ac(3; rUn - rrip) ~ (m„ - nip)—^ ~ 10 keV. (24) 

rriN 



This is in good agreement with the value from from Ref. 36|, Ac(3; mn—nip) = 14 keV. 



2. The electromagnetic interaction between nucleon magnetic moments, between the cur- 
rents associated with moving protons, and due to vacuum polarization. Those effects 
were calculated numerically to increase the mass difference by less than 30 keV (see 



e.g. 



561 ] where the various contributions where calculated with the AV18/IL2 model). 
Furthermore, the corrections to the Coulomb potential due to the proton's finite size 
have to be accounted for. In fact, these turn out to be the largest of the higher-order 
electromagnetic effects, since they modify the Coulomb energy by a fractional amount 
~ Tp/R"^, where R is the typical distance scale which dominates the Coulomb energy. 
Putting in r„ ^ 0.85 fm, R ^ 2.5 fm we might expect up to a 10% effect. Numerical 



evaluations [18[ however suggest a somewhat smaller number, reducing the mass dif- 
ference by about 33 keV. Importantly, the finite proton size decreases the impact of 
the electromagnetic interaction on the binding energy (thus increasing B(^He)) while 
the other effects listed above increase it (and so decrease B(^He)). Thus, each class 
of correction could individually induce an error of order 30 keV in our LO binding 
energy calculation, but they work in opposite directions, such that we can confidently 
say that their combined effect will not produce more than a 28 keV shift in A (3). 

Combining these two higher-order effects linearly with the 2 keV RGM uncertainty (see 
Appendix |A]) we find a potential higher-order correction due to long-distance effects which 
could be as large as Ac(3) = 44 keV. Employing Eq. ( 122|) . this produces an uncertainty of 

= 2.1 fm. (25) 
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The combined theoretical uncertainty in A (a) due to higher-order long- and short-range 
interactions is then, conservatively, taken to be 

^Lo = + <5,hon = 4.1 fm. (26) 

The central value for a„„ between the maximum adopted at A — > oo (see eq. ( 120|) ) and the 
minimal prediction at A = 800 MeV is —22.9 fm. EFT(/) thus yields, at leading order, a 
neutron-neutron scattering length of 

ann (EFT(^)) = -22.9 ± 4.1 fm. (27) 

Hence, the datum a„„ = —18.7 ± 0.7 fm is consistent with the input data set 
{onp,s5 c^^, B{d, t,^ He)} while the other datum a„„ = —16.1 ± 0.4 fm is inconsistent. 

C. The limit ann — ^ 

The interactions which generate a„„ and different: the short-range part is of the 

same structure but different in strength. In the nn case it is solely responsible for ann, 
whereas, for pp, it complements the Coulomb force to yield a^. It is therefore not obvious 
that for C"'" such that a„,„ ~ CLpp, the triton will be bound by the same amount as '^He is. 
Calculation at A(a) = (i.e. a„„ = a^) results in 

A(3) = -0.11 ±0.1 MeV, (28) 

i.e., two three-body systems of equal binding energy. The bound states result from the same 
short-range np interaction, but a purely short-range force between the like pair in the triton, 
and a combination of a similar short-range counter term plus the Coulomb interaction in 
Helium-3. The resultant approximate degeneracy in the binding energies is a reflection of 
the fairly small characteristic momenta. Comparing the nn and pp phase shifts resulting 
from interactions with A (a) ^ at A = 1.6 GeV (Fig. Hj), the nn interaction is found 
less repulsive than the pp one below E^^ ^ 0.65 MeV, equal around 0.65 MeV, and more 
repulsive for i^^.m. ~ 0.65 MeV. If bound states receive significant contributions from modes 
with k^^ ^ ^ymj^{2 MeV) =: Pbaiancc ~ 45 MeV, then the difference in Vnn and Vpp is 
naively expected to balance, yielding approximately the same binding energies, as we see in 
Eq. ( 128|1 . since ann ~ o-pp produces A(3) ~ 0. In fact, this analysis implies a momentum 
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distribution amongst the nucleons within the triton bound state which is dominated by 
momenta markedly smaller than the conventional estimate: Ptyp := -\/2 ■ 2/3mN ■ B{t) ^ 
100 MeV — at least as far as the momenta pertinent to the binding-energy difference A (3) 
are concerned. 

Furthermore, the reasoning implies that more deeply bound mirror nuclei will exhibit 
a larger difference in their binding energies, even though the respective two-body scatter- 
ing lengths are equal. In a world where not only A (a) = but additionally i?(^He) 3> 
i?(^He,exp), the uncharged mirror image would be not as deeply bound, B(t) <^ B{^B.e), 
as a result of the stronger repulsion of the uncharged 'neutrons' as compared to protons at 
relative momenta greater than about 20 MeV. 

This hypothesis is confirmed by the results of a RGM calculation (Fig. |5]). A smooth 
change of the three- nucleon parameter Asni increases -8(3) but leaves the two-nucleon sector 
invariant, i.e., ann = CLpp ~ —7.8 fm, Qnp^s ~ —23.75 fm, and the deuteron at its physical 
binding energy. The result is an increasingly less bound uncharged system relative to its 
charged mirror sibling, as conjectured above. In more detail, the triton is found not as deeply 
bound (intersection of dashed line with gray band in Fig. |5]) as in our world (black band in 
Fig. |5]) for a TNI producing the physical Helium-3 binding energy. This is a consequence of 
the more repulsive nn force at low momenta implied by A(a) = 0. Adjusting the TNI to 
yield a more deeply bound Helium-3 widens the gap between its ground state and the p — d 
breakup threshold because B{d) remains constant. The triton binding energy also increases 
in this procedure and is found larger than S(^He) for B{^B.e) ~ 12 MeV (A(3) > as shown 
by dashed line in Fig. |5]) becoming increasingly less bound relative to B{^B.e) as B{^B.e) 
increases further. This behavior of B(^B.e) is in accord with the qualitative discussion 
above. 

For systems with p^^^ much smaller than Pbaianco the opposite behavior is expected. If we 
assume, for the moment, that a„„ = a^, the above reasoning can be used to shed some 
light on the situation in the six-body system. In ^He, the two halo neutrons are very weakly 
bound together with an a core. Of course, it is possible that ^He would not be bound if 
ann were reduced to agree with a^. But, even if ^He were bound for this smaller a„„, the 
momenta in that bound state are low enough that the mirror nucleus ^Be with two 'halo' 
protons will not be bound. The pp interaction is more repulsive for those low momenta — due 
to Coulomb effects. An explicit calculation at an order in EFT(/) which generates ^He as 
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a shallow bound state is required to validate this hypothesis. At present, a leading-order 
analysis is not able to reach this level of accuracy. 

V. CONCLUSION 



The Coulomb energy (V^) of the three-nucleon system can be computed reliably with 



EFT(^) wave functions. In Ref. 30| an NLO EFT(/) computation gave a value of 660 ± 
30 keV for (Vc). In this work we have considered, in addition to Coulomb effects, the 
impact of the charge-symmetry-breaking NN operators which produce different (strong) pp 
and nn scattering lengths. We carried out the LO EFT(^) calculation for three different 
cutoffs A = 400, 800, and 1600 MeV using the modified renormalization-group method. We 
found a robust correlation between B{t) — 5(^He) and the difference of scattering lengths 



A(a) := ttr, 



"pp- 



The fact that this correlation is largely independent of the short-distance 



physics in the NN system indicates that three-nucleon operators do not contribute to this 
isospin-violating difference of binding energies at leading order. 

From the correlation and the experimental values of the tri-nucleon binding energy dif- 



ference and we infer: 



-22.9 ±4.1 fm. 



(29) 
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The uncertainty here has been assessed by adding hnearly estimates of the impact of higher- 
order, short-distance operators in the NN system and of neglected long-range effects {e.g., 
magnetic- moment interactions, as well as the nucleon mass difference). "Short-" and "long"- 
distance effects of higher order appear to contribute roughly equal amounts to the error bar. 

The result fl2U]) is due to operators that are first-order in isospin breaking, and so first- 
order perturbation theory with these operators, evaluated between charge-symmetric triton 
wave functions, could also have been employed (c.f. Ref. 4l|)- Here we performed an as- 
sessment of isospin violation in the Hamiltonian in which the relevant interactions were 
treated non-perturbatively. Such a calculation is more straightforward technically within 
the RGM. At the level of accuracy of our calculation, the only operator for which this non- 
perturbative/perturbative distinction might make a difference would be the pp Coulomb 
potential. But, even there, we anticipate that the second-order piece of the ^He Coulomb 
energy is of the same size as other effects neglected in this calculation. 

In order to refine the constraint ( !29|) it will be necessary to compute explicitly higher- 
order electromagnetic effects and the impact of the proton-neutron mass difference on the 
binding energies. Next-to-leading order and next-to-next-to-leading order triton wave func- 



tions in the charge-symmetric sector 30|, l32| should be considered. Analysis along the lines 
of Ref. 46| would also be needed so as to determine the order at which charge-symmetry- 



breaking three-nucleon operators enter the EFT(/) calculation. 
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Appendix A: Numerical stability 

The RGM is used to fit LECs and for predictions in the two- and three-body sector. In this 
section two analyses are presented to estimate the numerical uncertainty: the convergence 
of a B{t) calculation with respect to dimension and "quality" of the variational basis for 
a given set of LECs, and the dependence of C^^ on the parameters used to expand and 
regulate Coulomb functions in a Gaussian basis. 



23 



TABLE II: Uncertainty in the triton binding energy AB{t) due to the omission of specific ba- 
sis states, e.g., an addition of vectors of the deuteron- neutron grouping with particles in rela- 
tive S waves and support from approximately l/\/8 fm to infinity will change B{t) by less than 
0.0011 MeV. 



cluster 


[Sn <E) Spf'^ 


orbital angular momentum [li (X) ^2]^ 


Gaussian 


AB{t) [MeV] 


d-n 


Sl2 = l 


[0(5 


5 of 


71,2 e [0, 8] fm-2 


0.0011 


d-n 


Sl2=0 


[Od 


5 0]° 


71,2 e [0, 8] fm-2 


0.0003 








nf 




0.0003 






[2 a 


5 2]° 




0.0003 



In the first scenario, the interaction is specified through five LECs in Eq. (IHl): 
Ci,2 7 C*5"'^^, C*3Ni- Predictions for B(t) will depend on the dimension of the RGM basis, 
D = Di + D2 + -D3, and a "wise" choice of width parameters, = 1 . . . -Dj}, for each 

grouping i. We consider ourself wise because the widths are chosen to expand an object of 
limited size, which is estimated by the deuteron and triton binding energies to be described 
within a central potential whose range is set by the regulator cutoff A. A larger A relates to 
a shorter-range interaction, mandating larger widths '-fij to account for the larger values of 
the wave function resulting from the deeper well. Simultaneously, the exponential tail has 
to be modelled accurately by keeping the smaller widths corresponding to the longer-range 
part. In essence, larger As require larger bases but do not pose an in-principle limitation for 
the application of the RGM. To assess whether a certain variational basis expands the triton 
accurately, the supposedly complete basis, ,i = 1, . . . , Di} with D12 = 60 and D3 = 0, 
is extended by 5 vectors all taken from the dominant grouping — which, for the triton, is the 
deuteron-neutron one. In the new, {D' = {Di + 5) + -D2 + -D3)-dimensional space, the triton 
binding energy is calculated as a function of one width parameter, either 71 for the deuteron 
fragment, or 72 for the separation of the neutron from the deuteron. Figure [6] displays the 
graph of the function / (7j) = B{t, D') — B{t, D), with B{t, X) being the smallest eigenvalue 
of the system 

{(j)m\H\(f)n) = E{(f)m\(j)n) , (Al) 

with indices m,n specifying the X variational parameters Qj in Eq. ( ITSi) . The numerical 
uncertainties due to the finite basis are summarized in table [Tll In conclusion, a total 
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FIG. 6: Difference in B{t) due to the addition of basis vectors to a reference set - wl2wl2(0). This 
change is shown in the lowest panel as a function of either of the two width parameters 71.2 for a 
d-n (d-n) grouping in two S waves. The middle (top) panel displays the change for the nucleons 
in relative P (D) waves. 

uncertainty in the three-body binding energy -8(3) — the numbers are of the same order of 
magnitude for ^He — due to the truncation of the variational basis of 

A (RGM) = ±1.5 keV (A2) 

is assigned to this analysis. This value is markedly less than the naive 10 % LO EFT(/) of 
i?(t,exp), and hence the 120-dimensional basis is sufficient for the accuracy of this order. 

The LEG values depend on the loop regulator A and on the RGM basis that spans the 
space in which they are fitted to data. For a cutoff variation at A > 771^^ and a modification 
of the basis using states with support only for particle separations less than approximately 
l/m^r, both dependencies reflect a modification of non-observable high-energy modes ab- 
sorbed in the LEGs. However, the dependency on the model space is sought to be minimal, 
so as to allow application of and comparison with the LEG values found using other numer- 
ical methods. To assess this uncertainty, the 20-dimensional S-wave two-body basis used 
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throughout this work was refined in two ways: 

• add one basis state and determine C^^ as a function of the width parameter; 

• determine C^g as a function of the regulator parameter used for the irregular Coulomb 
function; 

In both cases, the change in the LEG was —0.01 < -jw- < 0.001, and the resulting effect on 
two- and three-body observables relevant for this work was small relative to the anticipated 
leading-order EFT(/) accuracy. 
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